Bose-Einstein supersolid phase for a novel type of momentum dependent interaction 
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A novel class of non-local interactions between bosons is found to favor a crystalline Bose-Einstein 
condensation ground state. By using both low energy effective field theory and variational wave- 
function method, we compare this state not only with the homogeneous superfluid, as has been 
done previously, but also with the normal (non-superfluid) crystalline phase and obtain the phase 
diagram. The key characters are: the interaction potential displays a negative minimum at finite 
momentum which determines the wavevector of this supersolid phase; and the wavelength corre- 
sponding to the momentum minimum needs to be greater than the mean inter-boson distance. 



Since Penrose and Onsager's first discussion [l| on ex- 
istence of "supersolid" , namely a phase with co-existence 
of superfluid and crystalline order, both experimental [2] 
and theoretical [3j attempts have been made for decades 
in the search of this novel phase. Recently reported ob- 
servation of "supersolid" phase in He-4 systems [J] re- 
vitalized this fundamental interest. Nevertheless, some 
subsequent experimental evidences as well as various pro- 
posed microscopic mechanisms [5] remain controversial. 

Progress on the physics of cold atoms and molecules 
opens a new possibility to study the "supersolid" phase 
thanks to clean and controlled experimental systems. 
One of the most fascinating facts is the unprecedented 
tunability of the interaction potentials due to internal 
degrees of freedom of atoms and molecules [fa, |7| , which 
allows one to address a theoretical question, namely what 
interaction potentials can support the supersolid phase in 
continuous space. Recent experimental progress on dipo- 
lar quantum gases allows to explore new physics of quan- 
tum many body systems with non-local interactions [8j- 
It is well established that non-local interaction potentials 
stabilize the supersolid phase on lattice [9| . The possibil- 
ity of finding a Bose-Einstein supersolid phase [lCf was 
also put forward for several continuum model systems 
such as dipolar quantum gases |ll| . atom-molecule mix- 
ture gases I12J1 an d Rydberg atom gases [7J, [l3| . Recently, 
Henkel et al [13J found that the Fourier transform of an 
isotropically repulsive van der Waals interaction poten- 
tial with a "softened" core has a partial attraction in 
momentum space, which gives rise to a transition from 
a homogeneous Bose-Einstein condensate to a supersolid 
phase due to roton instability. However, whether the su- 
persolid phase they found is stable against fluctuations 
and how it should compare with the non-superfluid (nor- 
mal) crystal phase has not been studied. Recent work 
on dipolar gases [14j showed that the dipolar dominating 
interaction does not support a supersolid phase in the 
phase diagram between the uniform superfluid and the 
normal crystal phase, where this phase had been specu- 
lated to exist. 

In this Letter, we show that interaction potentials, 
which display a minimum of negative value at a finite mo- 



mentum, lead to a modulating superfluid order, namely 
a Bose-Einstein supersolid (BES) phase. We perform ef- 
fective field theory analysis and variational calculation 
to determine not only the phase boundary between the 
uniform superfluid (USF) phase and BES, which has 
been previously analyzed by roton instability for dipo- 
lar [lfj, |ll| or van der Waals interaction [13[, but also 
the phase boundary between BES and the normal (non- 
superfluid) insulating crystal (IC) phase. We shall be- 
gin with a heuristic argument to show how a stripe BES 
phase should arise from the competition between kinetic 
and interaction energy in the regime of roton instability. 
Next, we shall study as a concrete example the "softened" 
dipolar interaction recently proposed for Rydberg atomic 
gases |3) . We will establish the ground state in the sense 
of variational principle and find a first order phase transi- 
tion from the uniform superfluid phase to the triangular 
crystalline BES phase. Finally, we shall compare the en- 
ergies of BES and IC phases of the same lattice configu- 
rations, and find a regime in which the triangular-lattice 
BES is stable and has lower energy than both USF and 
(normal) IC. The result is summarized in Fig. [TJ 

Hamiltonian and heuristic. To explore the physics of 
the BES phase, we start with the continuum Hamiltonian 
of two dimensional interacting bosons 
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where the first term of H corresponds to the kinetic en- 
ergy, and the second the two-body interaction energy. 

It is commonly accepted that the ground state for such 
a continuous bosonic system should be USF at the kinetic 
energy dominating regime. The USF phase is described 
by a coherent state |USF) = cxp(/ d 2 x^/ne i ' t ' ^ (x))\n) 
where n is the mean particle density, 4>o a constant phase, 
and |0) the vacuum state with no particle. The energy 
of this state is given by Stjsf = -jTlUQc = 0), where N 
is the mean particle number and i7(k) is Fourier trans- 
form of the interaction potential. We first analyze the 
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FIG. 1: The phase diagram of bosons with step-like interac- 
tion. Uniform superfhhd (USF), insulating crystal (IC), and 
Bose-Einstein supersolid (BES) phases are separated by 'solid 
lines' with "x" showing the data points from variational cal- 
culation. Analysis of the collective excitation spectrum shows 
the instability of USF at the 'red dotted line' and that of BES 
at the 'dark blue dash dotted line'. At low density, USF phase 
and IC phase exist; while at high density the new stable BES 
phase is found in the 'yellow shaded' regime. When n > 1, 
the IC state is not stable (see text). 



instability of the USF phase. This can be performed 
using an effective field theory approach [15|, [l6|. The 
real time action of this bosonic system is S[ip,?p] — 
J d 2 xdt{ih-ipdttp — Jtf'l'ipjip}}. Fluctuations on top of 
the uniform supcrfiuid state are considered by writing 
the boson field ip( x ,t) = [/°o + Sp^^e 1 ^, assuming Sp 
and |V<^>| are small. The quasiparticle spectrum is read- 
ily derived after integrating out the dp field: e(k) = 

V ^"(^T + 2n ^( k ))- For a potential that has a neg- 
ative minimum at a finite momentum, this spectrum at 
that momentum drops, eventually hits zero and becomes 
imaginary when increasing the density n. That suggests 
that the assumed USF (coherent) state is unstable to- 
wards possible crystalline order. 

To show the BES phase arises, we first give a heuris- 
tic argument by considering a simple stripe BES state 

|BES) = cxp (V3V(^&J, /2 + ^r&L Q/2 )) |fi>, where Q = 
[Q, 0], and Q the minimum point of E/(k). The energy of 
this state is given by E sBES = N {^- + \nU(Q)^ + 

E VSF . When the term 2ySL + |nJ7(Q) is negative, 
namely the interaction energy dominates over the kinetic 
energy, the stripe BES state has lower energy than the 
USF state. (We also go beyond the mean field state 
and compare with the two component fragmented state 
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are variational parameters |17j . and the coherent stripe 
BES state is found to have the lowest energy.) We thus 
conclude the BES state arises from the competition of 
kinetic energy and interaction energy. 
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FIG. 2: LEFT figure (a): Fourier transform of the step-like 
two-body interaction. LEFT figure (b): Shows the Bogoli- 
ubov quasiparticle spectrum for a USF state. The unit of 
£k is ft 2 /(mro). The plot shows the real part of the spec- 
trum with h — 1. The solid line corresponds Td — 5, the 
dashed line to r^ = 15, and the dashed dotted line r^ = 23. 
RIGBT figure: Shows the phase transition from the USF to 



the triangular crystalline BES phase. \4>k.\ 
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the occupation fraction of the lowest finite momentum. 



To be concrete, we further apply the two-particle inter- 
action of a step-like form V(r) = § if r < r^; V(r) = -§■ 

r .° . r 

otherwise. The form of this potential is an approximation 

to the interaction between polarized Rydberg atoms pro- 
posed in Ref. [7| . Two dimensionless parameters of this 
system are n = n x r 2 , and r^ = 2 "'^" — . n characterizes 
the relation between ro and the inter-particle distance, 
and Td characterizes the strength of interaction. A phase 
transition from USF to IC has been found when varying 
rd at the regime of h ss 0.9 [14j , The IC (single parti- 
cle per site) phase is described in a second quantization 
form by l^ic) = T[r c - |0)j where Ri is the direct lat- 
tice vector at site i, and the single particle wavefunction 
corresponding to c - is the Wannier function </>s (r). 

Ri tt* 

The Fourier transform of this step-like interaction is 
shown in FIG. 2(a). It is straightforward to obtain the 
excitation spectrum, which is shown in FIG. 2(b). It 
can be seen that the spectrum displays instability. The 
origin of this effect is that the Fourier transform of the 
interaction, C/(k), has a negative minimum at a finite 
momentum. Now the question is to find the stable varia- 
tional minimum in the coherent state space. With \G) = 

exp(/d 2 x(Hx)V>(x))|ft) (so that $(x)|G) = 0(x)|G», 
the energy of this state is readily given by: 
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where V(f) is the interaction potential. 

Variational analysis. We first check whether the sys- 
tem favors an extended or localized state. This purpose 
is fulfilled by applying the Gaussian ansatz which means 

cf)(r) = ^M-er^^ 1 . The total energy of this system is 
given by E t = E^ + Edi P , where the kinetic energy E^ = 



N 



2m<7' 



and the interaction energy E,u p = N 2 o^? ^3 (■?■)■ 



Here, g(x) is approximately 7r 2 (l 
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g(-). In thermodynamic limit 



N —$■ oo, interaction dominates and £ , (<t) = „ 2 :i g(^-). 

We found that as long as f <i 2 rV(r) > 0, er — > oo mini- 
mizes the energy, implying the system favors an extended 
state in space. Since for r^ > Jd 2 rV(f) > 0, we con- 
clude that the system favors an extended state when Vd 
is positive. 

Up to this point, what we know about this system 
is that it favors an extended state which is not neces- 
sarily uniform superfluidity. Having argued heuristically 
above that a momentum dependent interaction may fa- 
vor a finite momentum BEC, it is natural to compare 
the energy of a new, non-uniform coherent state which 
has discrete lattice symmetries. Thus, we can write the 
condensate wavefunction in such a form 4>{r) = (tp(r)) = 
V™J2k ^Ke' Kr with K = pG% + qG 2 , where Gi and 
G2 are two primitive vectors spanning the two dimen- 
sional reciprocal lattice. The corresponding ground state 
is \G) — exp (J2k V^V^k^k) 1^)- The or der parameter 
that characterizes the phase transition from uniform su- 
perfluidity to BE supersolidity is an occupation fraction 
at some finite momentum K, \4>k\ 2 = ^(^k&k)- 

In this assumed ground state subspace, the energy per 
particle is given by: 
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Now, the problem reduces to minimizing this energy func- 
tional with such a constraint £ K \4>k\ — 1; which is 
equivalent to the enforcement of conservation of the to- 
tal particle number. By 8 (E — p £ K ^k^k) / ' Wic) = ^' 
we obtain 
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where p is the chemical potential. 

We compute the energies for three different configura- 
tions — stripe, square and triangle lattices — and found 
that the triangular lattice is the most energetically fa- 
vored. The optimal lattice constant <zbes is found to 
be slightly larger than 2ir/Q m i n where U(Q m i n ) corre- 
sponds to the negative minimum of the potential. For 
the particular step-like interaction, Q min is related to r$ 
by Qmin ~ — ~ — • The transition between the uniform 
superfluid and the supersolid lattice is of first order as 
shown in FIG. H 

Dynamic stability of the Bose-Einstein supersolid 
phase. To study the stability of the BES phase, we 



further explore the fluctuations on top of the ground 
state. In the presence of BES phase, we can get the 
effective field theory for the density and phase fluc- 
tuations [l8j, 8p(x,t) and ip(jc, t), respectively, writing 

V>(x,i) = [p (x) + <5 /9 (x,i)] 1 /2 e ^(x,t) ; where po ( x ) = 

|</>(x)| 2 . The effective action for Sp and ip to quadratic 
order is 



S c s[Sp, if] — dt d xJzf(x, £) 



jsf = -h5pd^-h 2 (v-%) 2 -^p h 2 (Vv) 2 

8 ,/PO I 

d 2 xV(x-x / )<5p(x,t) ( 5 / 9 ( x / ,i), (5) 



where po( x ) is the modulus square of the condensate 
wavefunction. 

Since the effective theory possesses only discrete trans- 
lational symmetries, a Brillouin zone and its reciprocal 
lattice vectors can be defined. Let 77 (x) = P X- . The 
effective action in the momentum space is 
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+ c.c, (6) 
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C Kl K 2 (k) = ^^[(Ki+kJ^-Ka-kJlttKOKa-Kx-KfT) 
K 

where Ki,K2 are reciprocal lattice vectors and qk is 
the Fourier component of \J po (x) . The first Brillouin 
zone is divided into 'R' (right) and 'L' (left) subzones 
according to time reversal; the summation k £ R in 
Eq. ([B]) means summing over the 'R' subzone. The 
above effective theory is quadratic in fields. Formally, 
the action can be written in a block diagonal form as: 



Seff = £ k£R [ ^(k) pt(k)] 



B(k) T + B(ky 



T,(k) 



with 



A(k)*d t 
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where the crystal momentum k is a good quantum num- 
ber, and r?(k) and p(k) correspond to column vectors 
formed by {?7K+k} and {</9K+k} with K running over 
pGi + qG2- The energy spectrum is determined by the 
poles of S? k , i.e., det^ 1 ] =0 0. FIG.Hshows typical 




FIG. 3: The excitation spectrum on top of the BES state 
from the T(= (0,0)) point to X(= (^-£-,0)) with a BE s the 
lattice constant of the supersolid lattice. The energy unit is 
h /(mr ). There are branches of spectrum, the lowest three 
of which are shown above. 



results we obtained, which indicate the stability of the 
BES phase. However, for sufficiently large r<2, the spec- 
trum becomes imaginary near k = (crystal momentum) 
point indicating instability of the BES state. The stable 
regime of the BES state are shown in the phase diagram 

FiG.ru 

Variationally compared with insulating crystal. We 
further estimate the energy of the IC state. The IC state 
is described by I'J/ic) = Yin ^~ |0) where the single parti- 
cle wavefunction corresponding to the creation operator 
c- is the Wannier function <j)s(r). Here we consider 
the case where each localized wave function contains ex- 
actly one boson, forming a triangular lattice. The lat- 
tice constant a c is thus determined by the density ( a c — 
[2/(v / 3"-)] 1 ^ 2 ), which is different from the lattice constant 
of BES, <2beSi determined by the minimum point of f (k). 
The Wannier function is approximated by a localized 
Gaussian <j)^{r) ~ ~~r^ ex P(~ 2^ ) w ith aja c rang- 
ing from to 0.3 over which the overlap between neigh- 
boring Gaussian wavefunctions can be neglected [19j and 
the energy is obtained by Eic(<j) — (^icl-ff l^ic)- When 
n > 1, minimizing E\c(<j) gives a/a c — > 0.3, indicating 
the IC state with the given lattice constant is unstable. 

In the phase diagram (FIG. [T]), the BES is stable and 
is the most energy-favored in the 'yellow shaded' regime. 
The lower boundary is determined by comparing the en- 
ergy of BES state and Eic{<7 = 0.3a c ) while the right 
boundary is computed from the instability of the BES 
spectrum. In the 'unstable' regime, the proposed BES 
state has lower energy than the IC state but is not stable 
against quantum fluctuations. 

In conclusion, we studied a bosonic system with two- 
body interaction potentials which display a negative min- 
imum at a finite momentum and found a stable supersolid 
phase arising from Bose-Einstein condensation at finite 
momenta. The stability of this novel supersolid phase 
is checked against quantum fluctuation. A unique fea- 
ture of the BES state is that it breaks both U(l) and 
translational symmetry with a single order parameter, 
namely, the superfluid order parameter (ip (f)} is not 



only finite but also spatially modulated. The physical 
interpretation is that particles are not localized in space 
but condensed to a single, common wavefunction which 
is modulated like a solid. This is conceptually differ- 
ent from one of widely considered supersolid pictures of 
He-4 [lj, |j) in which supersolidity is a mixture of two or- 
derings: atoms form charge-density-wave order (a crystal 
structure) and in the same time vacancies or interstitials 
undergo usual (zero-momentum) Bose-Einstein conden- 
sation. For the conventional superfluid phase originated 
from zero-momentum BEC, there exists long range phase 
coherence but the phase correlation function is homo- 
geneous, not modulated, in space. In contrast, for the 
IC (insulating crystal) phase, particles are localized in 
space to each lattice site, so there is no long range phase 
coherence. Therefore, as prediction for cold gas experi- 
ments, a signature of the new BES phase is the modu- 
lated phase coherence. This new state also opens funda- 
mental questions for future studies, for example, how the 
supercurrent is affected by the simultaneous presence of 
crystalline ordering, and topological configurations such 
as a vortex coupled to a crystal defect. 

We thank Erhai Zhao and Uwe R. Fischer for helpful 
discussions. This work is supported by ARO Grant No. 
W911NF-07-1-0293. 

Note Added. Near the completion of this paper, a re- 
lated, independent study appeared [20J , which discovered 
by exact numerical algorithms for a similar model system 
a supersolid phase to occur in the same parameter regime. 
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